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Measurements of new physics at the forthcoming experimental program at CERN's Large 
Hadron Collider (LHC) will require a precise understanding of processes at next-to-leading order 
(NLO). This places increased demands for the computation of new one-loop amplitudes. This in 
turn has spurred recent developments towards improved calculational techniques. 

Direct calculations using Feynman diagrams are in general inefficient. Developments of more 
efficient techniques have usually centred around unitarity techniques where tree amplitudes are 
effectively "glued" together to form loops. The most straightforward application of this method, in 
which the cut loop momentum is in D = 4, allows for the computation of "cut-constructible" terms 
only, i.e. (poly)logarithmic containing terms and any related constants. QCD amplitudes contain, 
in addition to such terms, rational pieces which cannot be derived using such cuts. These "missing" 
rational parts can be extracted using cut loop momenta in D = 4 — 2e [pp. The greater difficulty of 
such calculations has restricted the application of this approach, although recent developments [||, 
0] have provided new promise for this technique. 

Recently the application of on-shell recursion relations ^ to obtaining the "missing" rational 
parts of one-loop processes ^\ has provided an alternative very promising solution to this prob- 
lem. In combination with unitarity methods an "on-shell bootstrap" approach provides an efficient 
technique for computing complete one-loop QCD amplitudes [^]. Additionally other new methods 
have also proved fruitful for calculating rational terms [§J. 

Such developments have again refocused attention on the optimisation of the derivation of 
the cut-constructable pieces of the amplitude. Deriving cut-constructible terms for any one-loop 
amplitude reduces to the computation of coefficients of a set of scalar bubble, scalar triangle and 
scalar box integral functions. Box coefficients may be found with very little work, directly from 
the quadruple cut of the relevant box function A unique box coefficient contributes to each 
distinct quadruple cut. Unfortunately triangle and bubble coefficients cannot be derived in quite so 
direct a manner. Multiple scalar integral coefficients appear inside a two-particle cut or triple cut. 
It is therefore necessary to disentangle the relevant bubble or triangle coefficients from any other 
coefficients sharing the same cut [jy, 10, 11]. The large number of NLO processes of interest for 



the LHC suggests that a completely automated computational procedure is highly desired. To this 
end we discuss, in this proceeding, a recently proposed method [[l2|, [l3|] for the direct, efficient and 
systematic extraction of bubble and triangle coefficients which is well suited to automation. 

1. Triangle coefficients 

Following in the spirit of the box coefficient [||] we would like to apply a triple cut to extract 
a triangle coefficient. Such a triple cut isolates a unique triangle coefficient but also contains 
contributions from scalar box integrals which share three of their four propagators with the triangle. 
The separation of the coefficient of a particular scalar triangle integral from any box coefficients 
can be effected by 



(1.1) 

f=0 



c = - [Inf f AiA 2 A 3 ] (t) 
Equation |0| instructs us to start by taking the triple cut of the desired triangle coefficient, 

ATA^A^) =A c ^ Cl+2 (-/,c 1 ,...,(c3-l),/i)A^ C3+2 H 1 ,C3,...,(c 2 -l),/ 2 ) 

xAf^H.a-.tci-l),!), (1.2) 
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menta. The cut momentum / of the triple cut depends on a single parameter t, and is parameterised 



in the specific form [ 1 1 



(/ '"'^ l + (7 2 -^ 2 ) ( ^ (/ l-( 7 2- 5l52 ),^ 1 + ^2 I- d-3) 

Here we have expressed the cut momentum in terms of a convenient basis of null vectors K[ and 

4 

jJw _ <-(5i/y)^ „b,ju _ ^-(5 2 /y)< 

Al " l-cwr 2 ) ' 2 " i-(^ 2 /r 2 ) ' ( j 

with two solutions for y = (K\ • K 2 ) ± VA with A = (Ki ■ K 2 ) 2 - K\K\. 

The Inf f is instructing us to series expand this cut integrand around t = oo. The t° component 
of this series expansion gives the desired triangle coefficient. For the three-mass case described 
above we must also average over the two solutions to y. Analytic continuation of / to complex 
momenta allows one- and two-mass triangles, containing three-point vertices, to be computed in a 
similar manner after setting the relevant masses in eq. (1.3) and eq. (1.4) to zero. In these cases 
only one solution to y survives. 

This procedure succeeds because of the specific momentum parametrisation we have chosen. 
The series expansion of the lnf t would in general give us rational coefficients a,- multiplying inte- 
grals over powers of t. Seen schematically this is 

a i j dtt'+ao Jdt + ai J dtt + ... + a max J dtt max , (1.5) 

and we would expect contributions to the scalar triangle coefficient from every term. It is easy to 
show though that all integrals over t will vanish, eliminating any such contributions. For example, 



jdtt~J Sl ^Jffi } ~ (4>-\£l\K\nVi + (K b 2 -\K 2 \K\--)^2 = 0, (1.6) 



,(K 2 -m 

with a similar result for other non-zero powers of t. The % are Passarino-Veltman reduction coef- 
ficients. 
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2. Bubble coefficients 

A similar procedure applies to the extraction of bubble coefficients. To compute the coefficient 
of a particular bubble we use a two-particle cut and must disentangle our desired coefficient from 
the scalar boxes and triangles which this will also contain. We start from the generic two-particle 
cut A?°A?°(t,y) = A^ Ci+2 (-/, ( Cl + 1), . . . ,c 2 , h ) ^(-/j, (c 2 + 1), . . . ,c h l) shown in fig- 
ure |l|(b), with l\ =l — K\. Having isolated a single bubble coefficient we parameterise the cut loop 
momentum /, which now depends upon two free parameters t and y, using 

(i-\=t(Kl>-\ + j(i-y)(x-\, = ~ << + \ + (x + V (2.i) 

This is expressed in terms of a basis of massless on-shell momenta K\ and %. # is an arbitrary 
free vector, which the final result is independent of, used to define k\'^ = — {S\ / y) with 

Y=(X ± \£i\X ± ) = (X ± W\\X ± ). 

The equivalent expression to eq. is then given by 

(2.2) 

»T(m) 



b = -i [Inf f pnfyMa] (y)] (f) t - \ £ [\nf t AiA 2 A 3 ] (t) 



The first term is the natural extension of the single Inf ? of eq. ( |1.1[ ) to the case of two free parame- 
ters. In this double series expansion we expand around y = °° and then t = 00 and again drop terms 
proportional to t, because the corresponding integrals disappear. Integrals over y are non- vanishing 
though and are related to the scalar bubble integral, Bq(K^), via / dyy m = Bq(K%) j(m + 1). 

The naively unexpected second term of eq. ( |2T2| ) involves a sum over all the triangles % that 
contain the original two-particle cut. Writing the two-particle cut integrand schematically in the 
form 

ao{t)+a\[t)y + ■ ■ . + a max {t )y + ^ - — — , (2.3) 

<^ tri hi \y yivi) 

allows us to understand why contributions from triangles arise. Solving (l(y,t)-K 2 ) 2 = 0, the 
additional propagator present in the triangles in the last term above, gives us yi(t). Inserting this 



solution into the momentum parameterisation of / given by eq. (gjj) leaves us with the momentum 



parametrisation of the triangles that we wish to separate from the bubble. This parametrisation 



differs importantly from eq. ( |1.3| ) in that the integrals over t do not vanish, as can be seen for 
example with 

Jdtt~J d 4 l iX ^ } ~ (x-\£i\K\n^ + (X-mK\'-)^2 * 0. (2.4) 

The remaining contributions to the bubble coefficient are then found by relating these non- vanishing 
integrals over t to scalar bubbles using 

n.) - (%)" •*>■- ( £y W (2 . 5) 



7 J A m \p t ""(^1-^2 

and T (0) = 0. The coefficients % n i are given by 
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